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ABSTRACT 

We study the clustering properties of the galaxy clusters detectable for the 
Planck satellite due to their thermal Sunyaev-Zel’dovich effect. We take the 
past light-cone effect and the redshift evolution of both the underlying dark 
matter correlation function and the cluster bias factor into account. A the¬ 
oretical mass-temperature relation allows us to convert the sensitivity limit 
of a catalogue into a minimum mass for the dark matter haloes hosting the 
clusters. We confirm that the correlation length is an increasing function of 
the sensitivity limits defining the survey. Using the expected characteristics of 
the Planck cluster catalogue, which will be a quite large and unbiased sample, 
we predict the two-point correlation function and power spectrum for differ¬ 
ent cosmological models. We show that the wide redshift distribution of the 
Planck survey, will allow to constrain the cluster clustering properties up to 
z « 1. The dependence of our results on the main cosmological parameters 
(the matter density parameter, the cosmological constant and the normalisa¬ 
tion of the density power-spectrum) is extensively discussed. We find that the 
future Planck clustering data place only mild constraints on the cosmological 
parameters, because the results depend on the physical characteristics of the 
intracluster medium, like the baryon fraction and the mass-temperature rela¬ 
tion. Once the cosmological model and the Hubble constant are determined, 
the clustering data will allow a determination of the baryon fraction with an 
accuracy of few per cent. 

Key words: galaxies: clusters: general - cosmology: theory - dark matter - 
large-scale structure of Universe - cosmic microwave background 


1 INTRODUCTION 


In the standard picture of structure formation based on 
the gravitational instability paradigm, clusters of galax¬ 
ies represent the largest gravitationally bound systems 
in the universe. This is the reason for their cosmological 
importance. In fact, since the expected displacements 
from their primordial positions are much smaller than 
their typical separations, clusters retain the imprint of 
the main cosmological parameters, which can hopefully 
be constrained by studying their properties. In the past 
this opportunity has been largely exploited using clus¬ 
ter counts and their redshift evolution (e.g. Eke et al. 
1998; Sadat, Blanchard & Oubkir 1998; Viana & Lid- 
dle 1999; Borgani et al. 2001; Verde, Haiman & Spergel 


2001), obtaining tight constraints on the matter density 
parameter flom and power-spectrum normalisation <rg. 

Galaxy clusters are also good tracers of the large- 
scale structure of the Universe. Their clustering signal 
can easily be detected, even with a relatively small num¬ 
ber of objects. In fact, they form in overdense regions of 
the cosmological density field and are strongly biased. A 
confirmation that the spatial distribution of galaxy clus¬ 
ters is highly correlated came from the statistical anal¬ 
ysis of the first extended optical surveys (Nichol et al. 
1992; Croft et al. 1997), which however may be affected 
by the presence of interloper galaxies. A better identifi¬ 
cation of galaxy clusters is possible in the X-ray band. 
Here their emission due to thermal bremsstrahlung from 
the hot intracluster plasma, is more centrally concen¬ 
trated because it depends on the square of the baryon 
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density. In the past years the completion of extended 
catalogues covering a large fraction of the sky allowed 
the computation of the clustering properties of X-ray se¬ 
lected clusters (Abadi, Lambas & Muriel 1998; Borgani, 
Plionis & Kolokotronis 1999; Moscardini et al. 2000a; 
Collins et al. 2001; Schuecker et al. 2001). The results 
obtained for different surveys are fairly different, reflect¬ 
ing a strong dependence on the cluster selection criteria 
(limiting flux or luminosity in a given band) adopted. 

It is important to notice that the clustering of 
galaxy clusters represents an important test for mod¬ 
els of structure formation. As shown by numerical sim¬ 
ulations, cluster evolution can easily be understood and 
interpreted, once they are identified as the most mas¬ 
sive dark haloes produced by gravitational collapse. As 
a consequence, one can model their clustering proper¬ 
ties starting from the dark halo properties, which can 
be obtained from an extension of the standard Press- 
Schechter formalism (e.g. Mo & White 1996). The re¬ 
sulting correlation function depends strongly on cosmol¬ 
ogy and the primordial power spectrum. In this paper 
we will use an improved theoretical model (Matarrese 
et al. 1997), paying particular attention to the redshift 
evolution of object clustering and to light-cone and se¬ 
lection effects. A treatment of the redshift-space dis¬ 
tortions is also included. This model has already been 
applied recently to compare the existing cluster data 
(both in the optical and in the X-ray bands) with model 
predictions, finding further evidence in favour of low- 
density models (Moscardini et al. 2000b; Moscardini, 
Matarrese & Mo 2001). 

In this context, a new window is now opening 
to extend the study of the cluster distribution in 
the microwave band, thanks to the thermal Sunyaev- 
Zel’dovich (SZ) effect, i.e. the change of the blackbody 
spectrum of the cosmic radiation background produced 
by Compton scattering off the hot intracluster medium. 
Observational data of some individual galaxy clusters 
have already been obtained, allowing a better morpho¬ 
logical study of individual objects. However, at present 
well-defined surveys required for statistical studies are 
not available yet, even though several ground-based ex¬ 
periments with this objective are being planned for the 
near future (e.g. BOLOCAM, AMIBA, CBI, AMI, etc.). 

The goal of a full-sky survey of SZ clusters will 
be reached only with space missions. In particular, the 
Planck satellite, which is scheduled for launch in 2007, 
will detect of order 10 4 clusters thanks to the optimal 
choice of filter bands around the frequencies relevant 
for the SZ effect. This large number of clusters, even 
if obtained with an angular resolution lower than that 
of the above-mentioned ground-based bolometer arrays 
and interferometers, will allow an accurate measurement 
of the cluster correlation function and constraints on its 
redshift evolution. 

This paper reports on an extended study of the 
clustering properties of galaxy clusters detectable for 
the Planck satellite. It is worth mentioning that an es¬ 
timate of the contribution from the correlation among 
clusters to the angular power spectrum of the cosmic 
microwave background radiation anisotropy due to fluc¬ 
tuations of the SZ effect has been obtained by Komatsu 


& Kitayama (1999) using an approach similar to the one 
applied here (see also Refregier et al. 2000). 

The plan of the paper is as follows. In Section 2 
we discuss the general characteristics of the objects de¬ 
tectable for Planck. Section 3 reviews the theoretical 
model for the correlation function of SZ galaxy clusters 
in the framework of different cosmological models. In 
Section 4 we discuss the implications of the previous re¬ 
sults and the possibility of constraining the main cosmo¬ 
logical parameters using the future Planck observations. 
Section 5 shows the dependence of the clustering predic¬ 
tions on the characteristics of the intracluster medium, 
like the baryon fraction and the mass-temperature rela¬ 
tion. Conclusions are drawn in Section 6. 


2 GALAXY CLUSTERS DETECTABLE 
FOR PLANCK THROUGH THERMAL 
SUNYAEV-ZEL’DOVICH EMISSION 


Galaxy clusters can be detected in the microwave regime 
due to their thermal Sunyaev-Zel’dovich effect. Hot 
electrons in the intracluster medium Compton-scatter 
the cold photons of the microwave background and re¬ 
distribute them towards higher frequencies. The result 
is a temperature decrement below 218 GHz, and an in¬ 
crement above. 

The SZ effect in direction 9 is quantified by the 
Compton-?/ parameter, 

_ kn f 

y{9) =-% / dl n e (0,l) T(0,l) , ( 1 ) 

meC- 1 J 


where T is the electron temperature, n e the three- 
dimensional thermal electron density (hereafter as¬ 
sumed to follow a King profile), and a T is the Thomson 
scattering cross section. 

Assuming an isothermal distribution of gas and ne¬ 
glecting the background noise (see the following discus¬ 
sion), the total SZ signal received from a cluster is sim¬ 
ply the integral over the solid angle covered, i.e. 


(2) 

where Dd is the angular-diameter distance to the clus¬ 
ter, which of course depends on cosmology. The total 
number of electrons in the cluster N e can be assumed 
to be proportional to the virial mass M, 


r _ 1 + fa . M 

e ~ 2 Jb m ’ 


where m p is the proton mass, fb is the baryon fraction of 
the cluster mass and }h is the hydrogen fraction of the 
baryonic mass (here assumed as 76 per cent). Assuming 
an isothermal gas distribution in virial equilibrium, it 
is possible to relate the cluster temperature (in keV) to 
its mass, 


T = 6.03 


M 

ItruFi Mq 


2/3 

E 2/3 (z) 


Avir(^-) 

178 


( 4 ) 


The quantity A v i r (z) is the mean density (in units of 
the critical density at redshift z) of the virialised halo 
computed from the spherical collapse model; E(z) is 
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the Hubble constant at redshift z in units of its present 
value, E(z ) = H(z)/Ho. The normalisation of the mass- 
temperature relation is taken from the analysis of hydro- 
dynamical simulations by Mathiesen & Evrard (2001). 
Using this set of relations, the quantity Y depends on 
the cosmological model, and on cluster mass and red- 
shift only, i.e. Y = Y(M,z). 

For modelling the cluster population detectable for 
Planck we have to take into account the characteristics 
of the satellite detector (e.g. Haehnelt 1997; Bartlett 
2000; Bartelmann 2001). Its relatively poor angular res¬ 
olution will not allow most of the cluster population 
to be resolved (e.g. Aghanim et al. 1997; Hobson et al. 
1998). Unresolved clusters below the detection limit will 
produce a Compton-?/ background 7/bg. Neglecting back¬ 
ground correlations, the average background fluctuation 
level can be estimated as 


( A l/bg ) 2 


dz 


dV_ 

dz 


U +*) 3 


dMn(M,z) Y 2 (M,z)( 5) 


(e.g. Bartelmann 2001). In the previous expression, dV 
is the cosmic volume per unit redshift and unit solid 
angle, and n(M, z) is the cluster mass function which 
can be estimated using the Press-Schechter approach or 
more recent extensions, see below. 

Moreover, the original cluster SZ signal will be con¬ 
volved by the Planck beam profile w(6), which depends 
on the observing frequency. We assume for our purposes 
that cluster selection will be mainly performed at the 
highest frequencies and approximate the beam profile 
with a Gaussian with r.m.s. a w = 5 arcmin. 

We assume that a cluster is detectable for Planck 
if its integrated, beam-convolved Compton-?/ parameter 
is larger than a given sensitivity limit, i.e. Y > U m i n , 
where 


Y = 



d 2 e 


d 2 d' y{9') exp 


(9-9) 2 
2 ct 2 


( 6 ) 


The external integral covers the solid angle where the 
signal is large enough to be detected. Using its nominal 
temperature sensitivity and the solid angle of the beam, 
the sensitivity limit for Planck can conservatively be set 
to V'min = 3 x 10 -4 arcmin 2 (Bartelmann 2001). 


3 THE CLUSTERING MODEL 

Our model for predicting the clustering properties of the 
SZ galaxy clusters derives from a method already ap¬ 
plied to the study of clusters detected in the X-ray and 
optical bands (Moscardini et al. 2000a,b; Moscardini, 
Matarrese & Mo 2001). We will give a brief description 
of the technique only and refer to the original papers 
for a more detailed discussion. 

Matarrese et al. (1997; see also Moscardini et al. 
1998; Yamamoto & Suto 1999; Suto et al. 2000; Hamana 
et al. 2001) developed an algorithm for describing the 
clustering on our past light-cone taking into account 
both the non-linear dynamics of the dark matter dis¬ 
tribution and the redshift evolution of the bias factor. 
The final expression for the observed spatial correlation 
function £ 0 b s in a given redshift interval Z is 



Y . (arcmin 2 ) 

min ^ ' 


Figure 1. The dependence of the correlation length ro 
on the sensitivity limit y m i n (in arcmin 2 ) is shown for 
AC DM (solid lines), OCDM (dotted lines) and rCDM mod¬ 
els (dashed lines). Heavy lower and light upper lines refer 
cluster samples with 2 : < 0.3 and 2 : > 0.3, respectively. 


£obs if ) 


f z dz 1 dz2J7(zi)JJ'(z2) £obj(nzi,Z2) 
[f z dz 1 JI{z 1 )] 2 


(7) 


where J\f(z) = A f{z)/r(z), Af{z) is the actual redshift 
distribution of the catalogue and r(z) describes the re¬ 
lation between the comoving radial coordinate and the 
redshift. In (li]), £obj {r, zi, Z2) represents the correlation 
function of pairs of objects at redshifts zi and 22 with 
comoving separation r, which can be safely approxi¬ 
mated as £obj(r, 21, 22) « b e ff (zi)b e s ( 2 2 )£ m (r, 2 aV e). Here 
£ m is the dark matter covariance function and 2 aV e is a 
suitably defined intermediate redshift (see Porciani 1997 
for a discussion of possible choices). 

A fundamental role in the previous equation is 
played by the effective bias 6 e ff, which can be expressed 
as a weighted average of the ‘monochromatic’ bias fac¬ 
tor b(M, z) of objects with some given intrinsic property 
M (like mass, luminosity, etc): 


b ea (z) = JV(z)- 1 dlnM' b{M',z) N(z,M'), (8) 

J M 

where A f(z, M) is the number of objects actually present 
in the catalogue with redshift within dz of 2 and mass 
within dM of M, whose integral over In M is A fz). 

In the case of galaxy clusters, it is possible to fully 
characterise their properties by the mass M of their 
hosting dark matter haloes at each redshift 2 . This is 
consistent with the hierarchical model of structure for¬ 
mation where clusters are expected to form by merg¬ 
ing of smaller mass units. The comoving mass function 
n(z, M) of dark matter haloes, required both in eq. (|j) 
and in the computation of J\f(z,A4), can be estimated 
using some extension of the Press-Schechter formalism. 
Moreover, we can adopt the Mo & White (1996) relation 
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(or improvements thereof; see Sheth & Tormen 1999) for 
the monochromatic bias required in eq. (|[). More pre¬ 
cisely, we adopt the analytic relations obtained by Sheth 
& Tormen (1999; see also Sheth, Mo & Tormen 2001) 
which are very well reproduced by the results of high- 
resolution A-body simulations (Jenkins et al. 2001). 

For incorporating the fully non-linear regime, we 
use the fitting formula by Peacock & Dodds (1996) 
which allows the analytic computation of the redshift 
evolution of the dark matter covariance function £ m . 
We also include the effects of redshift-space distortions 
using linear theory and the distant-observer approxima¬ 
tion (Kaiser 1987). 

Finally, in order to predict the abundance and clus¬ 
tering of galaxy clusters detected through their thermal 
SZ emission, we need to relate the expected character¬ 
istics of the samples (which in the case of SZ detected 
objects can be expressed in terms of the sensitivity limit 
Vmin) to a corresponding halo mass at each redshift. 
Such a relation is provided by the mass-temperature re¬ 
lation (eq. [ij). 

Changing the sensitivity limits, the characteristics 
of the detected cluster population change. In deep cata¬ 
logues, faint objects (corresponding to low-mass haloes) 
can be detected. Due to the strong mass dependence 
of the monochromatic bias, a different clustering am¬ 
plitude is expected. In Fig. |lj we show the dependence 
of the correlation length ro (defined as the scale where 
Cobs is unity) on the limits used to define the SZ sur¬ 
veys, i.e. the limiting sensitivity T m i n . The results are 
shown for galaxy clusters having z < 0.3 (heavy lower 
lines) and z > 0.3 (light upper lines) because at this 
redshift the expected Planck sample is roughly divided 
into two subsamples with a similar number of objects. 
We consider three different cosmological models, whose 
parameters are listed in Table |I| We find that ro is an 
increasing function of the limits defining a cluster sur¬ 
vey: the larger T m in is, the more clustered the objects 
are, independent of the cosmological model. It is inter¬ 
esting to note that a similar trend was also predicted 
for galaxy clusters detected in the X-ray (Moscardini et 
al. 2000b) and optical bands (Moscardini et al. 2001). 


4 COSMOLOGICAL IMPLICATIONS 

We now discuss how the clustering properties of the SZ 
clusters depend on the cosmological parameters. We set 
the SZ sensitivity limit to that expected for the Planck 
satellite, i.e. V m i n = 3 x 10 -4 arcmin 2 . We will consider 
in detail three different cosmological models. They have 
a cold dark matter (CDM) power spectrum, a primor¬ 
dial power spectral index n = 1, and a shape param¬ 
eter r = 0.21 in agreement with an extended set of 
observational data (e.g. Peacock & Dodds 1996). More¬ 
over, they have the power-spectrum normalisation ob¬ 
tained by Viana & Liddle (1999) to reproduce the local 
cluster abundance. Notice that this normalisation has 
been challenged by very recent analyses, which, even 
if based on quite different approaches, seem to con¬ 
verge to substantially smaller values for erg (Reiprich 
& Bohringer 2002; Viana, Nichol & Liddle 2002; Sel- 


Table 1. The parameters of the cosmological models. Col¬ 
umn 2: the present matter density parameter Pom; Column 
3: the present cosmological constant contribution to the den¬ 
sity Poa; Column 4: the primordial spectral index n; Column 
5: the local Hubble parameter h ; Column 6: the shape param¬ 
eter F; Column 7: the spectrum normalisation erg-, Column 
8: the baryon fraction fg. 


Model 

^Om 

UoA 

n 

h 

r 

<78 

fb 

ACDM 

0.3 

0.7 

1.0 

0.7 

0.21 

0.99 

0.128 

OCDM 

0.3 

0.0 

1.0 

0.7 

0.21 

0.84 

0.128 

rCDM 

1.0 

0.0 

1.0 

0.5 

0.21 

0.56 

0.212 


jak 2001). Combined analyses of data coming from the 
2dF Galaxy Redshift Survey and from measurements 
of the cosmic microwave background anisotropies seem 
to corroborate these low values of the power spectrum 
normalisation (Lahav et al. 2002). We will discuss later 
the effect of changing the assumed value for erg. The 
models correspond to different geometries. Specifically, 
we consider a flat model with matter density param¬ 
eter Ho m = 0.3 and cosmological constant (hereafter 
ACDM); an open model with flom = 0.3 and vanishing 
cosmological constant (OCDM); and an Einstein-de Sit¬ 
ter model (rCDM). The two low-density models have a 
local Hubble constant h = i7o/(100km/s/Mpc) = 0.7, 
while we assume h = 0.5 for the rCDM model. The 
baryon fraction is set to fb = 0.075h~ 3 / 2 , as suggested 
by the analysis of X-ray data made by Mohr, Math- 
iesen & Evrard (1999). Table [I] summarises the model 
parameters. 

In the left panel of Figure we show how the prop¬ 
erties of the cluster population detectable for Planck 
change as a function of redshift. Locally objects with a 
relatively small mass will be detectable, of the order of 
galaxy groups: in fact M m i n Js 10 14 h _1 Mq for z <, 0.1, 
independent of the cosmological model. At higher red- 
shifts, the minimum mass increases: at z ~ 0.5 the 
Planck catalogue will include clusters with mass larger 
than M m i n ~ 10 14 ' 6 h -1 Mq for low-density models and 
M m [ n « 10 14 ' 3 h _1 Mq for the Einstein-de Sitter model. 
Note that the differences between ACDM and OCDM 
are almost negligible. In fact the SZ clusters detected by 
Planck have relatively small redshifts, consequently the 
presence of the cosmological constant has small impact 
on the distance relations and the mass-temperature re¬ 
lation. 

The right panel of Figure ^ presents the redshift de¬ 
pendence of the effective bias factor 6 e g. As expected, 
the bias is a strongly increasing function of z. This is 
due to two different reasons. First, the bias is itself an 
increasing function of mass, and we found an increase of 
Afmin with redshift; second, even fixing the halo mass, 
the (squared) bias is inversely proportional to the clus¬ 
tering of the underlying dark matter distribution which 
is growing with time. This also explains the largest val¬ 
ues of 6 e fi obtained for the rCDM model: in fact the 
growth factor is larger in an Einstein-de Sitter cosmol¬ 
ogy than in low-density models. 
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Figure 2. The minimum mass (left panel) and the effective bias factor (right panel) for clusters detectable for Planck are 
plotted as functions of redshift. Different line types refer to different cosmological models: A CDM (solid lines), OCDM (dotted 
lines), tCDM (dashed lines). The model parameters are summarised in Table [1 


Figure presents the predicted two-point correla¬ 
tion function £ 0 b s computed for cluster samples hav¬ 
ing z < 0.3 and z > 0.3 (lower and upper lines, re¬ 
spectively). In the plot ‘mock-observational’ 1 -a error 
bars are given (only for the ACDM model, for clarity). 
They are obtained by bootstrap resampling the num¬ 
ber of expected pairs in each separation bin (Mo, Jing 
& Borner 1992). £ 0 b s for the OCDM model is almost 
indistinguishable from the A CDM case. This confirms 
that the clustering of SZ clusters is almost insensitive 
to the cosmological constant, as already found for ob¬ 
jects detected in the X-ray and optical bands (see also 
the discussion of the following Figure^). The predicted 
correlation function for the rCDM model is lower by 
approximately a factor of two, both for z < 0.3 and 
z > 0.3. The correlation lengths ro (still defined as the 
scales where £ 0 b s is unity) for the low-redshift catalogues 
are 15.9 ± 0.4, 15.7± 0.6 and 11.9 ± 0.2 h -1 Mpc for the 
ACDM, OCDM and rCDM models, respectively. As ex¬ 
pected from the behaviour of the bias factor displayed 
in Figure we find that the clusters are more clus¬ 
tered at higher redshift: the corresponding values of ro 
for the sample with 2 > 0.3 are 22.8 ± 0.7, 22.7 ± 1.2 
and 17.3±1.1 h -1 Mpc, for ACDM, OCDM and rCDM 
models, respectively. Again, 1 -a bootstrapping error es¬ 
timates are given. 

The differences expected between low- and high- 
density models are confirmed by the power-spectrum 
Pobsijt ), which is shown in Fig. fu. Results are given for 
the cluster sample with z < 0.3 only. Again the dif¬ 
ferences between open and flat models with fiom = 0.3 
are negligible, while the power-spectrum for the rCDM 
model is approximately a factor of two lower. 

In order to better quantify the redshift evolution of 
clustering, we show in Fig. ||the values of the correlation 
length computed using redshift bins of size Az = 0.2 



r (hr 1 Mpc) 


Figure 3. The spatial correlation function £ 0 bs( r ) is shown 
for clusters detectable for Planck with 2 : < 0.3 (lower heavy 
lines) and z > 0.3 (upper light lines). The results refer to 
ACDM and tCDM models (solid and dashed lines, respec¬ 
tively). Results for OCDM are not shown because they are 
almost indistinguishable from the A CDM model. Error bars 
(shown only for ACDM ) are 1-cr bootstrap estimates. 


for the three cosmological models. From z ~ 0 to z ~ 
1, ro changes by a factor of almost two. Locally the 
SZ clusters detected by Planck are expected to have 
a correlation function only slightly larger than (local) 
galaxy groups or poor (optical) clusters. The predicted 
clustering signal at high redshift is quite large (albeit 


© 0000 RAS, MNRAS 000, 000-000 











6 Moscardini et al. 




Figure 4. Theoretical predictions for the power-spectrum 
P a bs(^) are displayed for the Planck sample of SZ galaxy 
clusters. Results are presented for KCDM (solid line), 
OCDM (dotted line) and tCDM models (dashed line) and 
refer to clusters with z < 0.3 only. 


more uncertain due to the smaller number of objects) 
and can certainly be measured, at least in the case of 
low-density models. Note that, given its very large area, 
the Planck survey will be one of the best opportunities 
to measure the redshift evolution of cluster clustering. 
An alternative possibility (but in the X-ray band) will 
be the survey proposed for the XMM/Newton satellite 
(Refregier, Valtchanov & Pierre 2001). In this case the 
expected limiting flux will be approximately Sii m = 5 x 
10- 15 erg cm -2 s _1 in the 0.5-2 keV band, but over a 
limited area of 64 square degrees. 

So far, we considered models with fixed cosmologi¬ 
cal parameters. Now it is interesting to discuss the pos¬ 
sibility of constraining these parameters using the clus¬ 
tering properties of Planck clusters. We start by con¬ 
sidering the changes of the correlation length ro when 
the matter density parameter flom is varied. The re¬ 
sults are presented in Figure for galaxy clusters with 
2 < 0.3 and 2 > 0.3 (heavy and light lines, respectively). 
All models shown here have the spectrum normalisation 
required to reproduce the local cluster abundance. We 
use again the relation found by Viana & Liddle (1999), 
namely os = O.56S2(0)(", where C = 0.34 and C = 0.47 
for open and flat models, respectively. Different choices 
for the normalisation will be discussed later. Primordial 
spectral index, shape parameter, Hubble parameter and 
baryon fraction are taken from Tab. [jj. We find that 
varying flom from 0.2 to 1.0 changes ro by 50 per cent 
for clusters with 2 < 0.3, and 30 per cent for clusters 
with 2 > 0.3. More precisely, the predicted correlation 
length decreases as fiom increases. Moreover, we notice 
that the cosmological constant increases the correlation 
length by less than 5%. In particular, at high redshifts 
the curves for flat and open models are almost iden¬ 


Figure 5. The redshift evolution of the correlation length ro 
is shown for AC DM (open circles), OCDM (open squares) 
and tCDM models (filled squares). Results are obtained di¬ 
viding the expected cluster samples in redshift bins with size 
Az = 0.2. Error bars (for clarity only shown for the A CDM 
and tCDM models) are 1 -o bootstrap estimates. 


tical. This is due to the relatively large value of the 
minimum mass of the detected clusters which implies 
the inclusion of objects with very similar properties. A 
similar result, i.e. the impossibility of constraining the 
presence of the cosmological constant using the cluster¬ 
ing of galaxy clusters, was also found analysing cluster 
data in the optical and X-ray bands (Moscardini et al. 
2000b; Moscardini, Matarrese & Mo 2001). 

In Figure ^ we show how the correlation length ro 
changes when the spectrum normalisation <rg is varied. 
Again, results are shown for galaxy clusters with 2 < 0.3 
and 2 > 0.3 (heavy and light lines, respectively). In the 
low-redshift bin, ro shows a small dependence on os and 
no significant differences are found using the Viana & 
Liddle (1999) normalisation compared to the more re¬ 
cent values obtained by Reiprich & Bohringer (2002), 
Seljak (2001) and Viana, Nichol & Liddle (2002). How¬ 
ever, for clusters with 2 > 0.3 the variation is expected 
to be quite large. The predicted correlation length for 
low-density models changes from ro « 35 to ro ~ 20h -1 
Mpc as os is increased from 0.3 to 1.5. For the rCDM 
model the change is more limited: ro « 23 h~ x Mpc for 
og = 0.3 and ro ~ 14 — 15/V 1 Mpc for 0.7 <rg < 1.5. 
In particular, when the normalisation is changed from 
the largest value (Viana & Liddle 1999) to the lowest 
one (Viana et al. 2002), ro increases by roughly 20 per 
cent. 

We also checked the effect on cluster clustering 
when only the Hubble parameter h is varied, keeping 
all other parameters fixed as in Tab. [I] (the results are 
not shown in the Figure). Varying h from 0.4 to 1, the 
correlation length changes at most by 10%, independent 
of the cosmological model and redshift bin. This varia- 
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Figure 6. The predicted correlation length ro is plotted as 
a function of the present matter density parameter fiom- Re¬ 
sults are presented for CDM models with T = 0.21 and ag 
chosen to reproduce the local cluster abundance (Viana & 
Liddle 1999). Solid lines: flat cosmological models (i.e. with 
non-zero cosmological constant); dotted lines: open models 
with vanishing 1 2 qa . Heavy lower and light upper lines refer 
to galaxy clusters with z < 0.3 and z > 0.3, respectively. 

tion is certainly too small to allow any constraints on 
the Hubble constant. 


5 DEPENDENCE ON THE PROPERTIES 
OF THE INTRACLUSTER MEDIUM 

The SZ signal of galaxy clusters strongly depends on 
the properties of the intracluster medium. In fact the SZ 
signal reflects the distribution of hot gas in the clusters, 
as shown by eq. (ill). So far, we employed a simple model: 
a baryon fraction fixed by the X-ray data of Mohr et 
al. (1999), a King profile for the gas distribution and 
a virial scaling relation between mass and temperature. 
In the following subsections we will study how sensitive 
the previous results are to these assumptions. Moreover, 
we will discuss if the clustering of SZ galaxy clusters 
can be used to directly constrain the properties of the 
intracluster medium. 

5.1 Baryon fraction 

In Figure ^ we show the dependence of the predicted cor¬ 
relation length ro on the baryon fraction fb- We consider 
cluster samples with z < 0.3 and z > 0.3 (upper and 
lower panels, respectively), and two different cosmolog¬ 
ical models: A CDM and tCDM (left and right panels, 
respectively). We do not present results for OCDM be¬ 
cause of their similarity with the KCDM model. Three 
different relations between fb and the present Hubble 
parameter h are assumed. First, we assume a perfect 
knowledge of the value of the Hubble constant: h = 0.7 



Figure 7. The predicted correlation length ro is shown as a 
function of the spectrum normalisation ag. Heavy lower and 
light upper lines refer to clusters with z < 0.3 and z > 0.3, 
respectively. Results are presented for KCDM (solid lines), 
OCDM (dotted lines) and tCDM models (dashed lines). Val¬ 
ues of ro correspond to different estimates of ag obtained by 
Viana & Liddle (1999), Seljak (2001), Reiprich Sz Bohringer 
(2002) and Viana, Nichol & Liddle (2002) are shown by solid 
circles, solid squares, open triangles and open squares, re¬ 
spectively. 

and h = 0.5 for KCDM and tCDM, respectively. Even 
if the second value appears to be too low when com¬ 
pared to recent observational estimates (e.g. Freedman 
et al. 2001), it is required in the case of an Einstein- 
de Sitter model to avoid a strong conflict between the 
age of the Universe and that of globular clusters. Sec¬ 
ond, we vary h in agreement with the constraints from 
primordial nucleosynthesis: fb = 0.024/flomh 2 (see e.g. 
the review by Schramm 1998). Finally we adopt again 
the relation given by Mohr et al. (1999) from the analy¬ 
sis of the X-ray data of 45 galaxy clusters. Notice that, 
given fb, the values of h corresponding to the previous 
relations can vary quite substantially. Consequently, the 
predicted scatter of ro indicates how well fb can be con¬ 
strained with these clustering data. 

In general we find small changes when we adopt 
the relations from primordial nucleosynthesis and X-ray 
data: the predicted values of ro differ only by 0.5-1 h _1 
Mpc for 0.05 £ fb m 0.25, independent of the cosmolog¬ 
ical model. Thus, even with a well-determined correla¬ 
tion length, as expected with the Planck survey, only a 
rough estimate of the baryon fraction will be possible. 
The situation changes if we assume independent knowl¬ 
edge of the Hubble parameter h. Considering the ex¬ 
pected l-er error bars on ro for the low-redsliift dataset, 
fb can be estimated from clustering data with typi¬ 
cal uncertainties smaller than 0.02-0.03 for the KCDM 
model and 0.05 for the tCDM model. At higher red- 
shifts (z > 0.3), one would expect a stronger dependence 
of the correlation length on the baryon fraction. In fact, 
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Figure 8. The predicted correlation length ro of galaxy clusters with z < 0.3 (upper panels) and z > 0.3 (lower panels) is 
plotted as a function of the baryon fraction fb- Results are presented for the ACDM (left panels) and rCDM models (right 
panels). Different lines refer to different assumptions on the relation between the baryon fraction and the Hubble parameter h: 
constant h, equal to 0.7 and 0.5 for ACDM and tCDM models, respectively (solid lines); /{, = 0.024/Tio ni /i 2 as predicted from 
primordial nucleosynthesis (dotted lines); fb = 0.075h —3 / 2 as derived from X-ray data (dashed lines). The dashed regions show 
the 1 -a ranges for ro as obtained using the parameters from Tab. 111. 


if the value of fb is increased, clusters with lower mass 
are included in the Planck catalogue. Since the mass 
function is so steep, almost all clusters will have masses 
near the lower mass limit, which essentially determines 
the bias factor and consequently the correlation length. 
This is confirmed by our results: the change of ro with fb 
at high redshift is quite substantial (about 40 per cent 
across the plots). However, the expected 1-u error bars 
in the determinations of the correlation length increases 
in this redshift bin and the possibility to constrain fb 
remains roughly similar to the z < 0.3 dataset. 


5.2 Dependence on the parameters of the 
mass-temperature relation 

The results obtained so far were based on a fixed mass- 
temperature relation (eq. ji|). The underlying assump¬ 
tion is an isothermal distribution of gas in virial equi¬ 
librium, which is in discrete agreement both with nu¬ 
merical simulations (e.g. Bryan & Norman 1998; Frenk 
et al. 1999; Mathiesen & Evrard 2001) and observational 
estimates (Xu, Jin & Wu 2001; Finoguenov, Reiprich & 
Bohringer 2002). However the parameters in the relation 
are not well known. For example, different groups ob¬ 
tained different numerical values for the proportionality 
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constant. Moreover, the presence of radiative cooling in 
cluster simulations changes both the normalisation and 
the slope of the mass-temperature relation (e.g. Muan- 
wong et al. 2001). Finally also the redshift evolution 
of the relation, which cannot be constrained by present 
observational data, can be modified by different mech¬ 
anisms, such as for instance the presence of an entropy 
floor (e.g. Tozzi & Norman 2001). 

For testing the robustness of our previous results we 
will now study the variation of the correlation length tq 
in response to changes of the parameters in the mass- 
temperature relation. We modify eq. (Q) as follows: 


T = a 1 


M 


10 15 /i _1 Mo 


' {l+z) a 3E 2/3 (z) 


A vi r (z j 


178 


1/3 


.(9) 


Here the parameters a 1: a 2 and a 3 represent the nor¬ 
malisation, the slope of the relation and the exponent 
of the redshift dependence, respectively. The results 
of the previous section correspond to (a 1 ,a 2 ,a 3 ) = 
(6.03,2/3,0). The cosmological model is fixed to the 
A.CDM model (see Table [ij). 

Figure shows the variation of the correlation 
length when only one parameter is changed. Results for 
clusters with a < 0.3 and a > 0.3 are presented in the 
left and right panels, respectively. For showing all curves 
in the same plots, we normalise the abscissa range us¬ 
ing a parameter A defined as A = (a; — a*)/(a/ — a\), 
where a\ and of are the lower and upper limits of 
the considered range for the parameter cu. We consider 
the following intervals which are expected to cover the 
whole parameter region suggested by theoretical and 
numerical models: 3 < aq/keV < 10, 0.3 < a 2 < 0.9, 
— 1 < a 3 < 1. As expected, the redshift dependence (a 3 ) 
is absolutely negligible in the low-redshift sample, but 
it also remains very weak for clusters with z > 0.3. The 
results are different when changing a 1 and a 2 . Increas¬ 
ing the normalisation of the mass-temperature relation 
from 3 to 10 keV, ro is decreased from approximately 
18 to 14.5h” 1 Mpc for z < 0.3 and from approximately 
26 to 21 h -1 Mpc for a > 0.3. Similarly, when the slope 
a 2 changes from 0.3 to 0.9, the correlation length in¬ 
creases from « 13 to « 17 h~ l Mpc for 2 < 0.3 and 
from as 21 to ss 23.5/i -1 Mpc for z > 0.3. We recall 
that the predicted 1-a error bars for r o are 0.4 and 0.7 
for the cluster samples at low- and high-redshift. Thus, 
using the clustering properties of SZ galaxy clusters to 
constrain cosmological parameters is hampered by our 
ignorance of the mass-temperature relation. 

However, the dependence of ro on the details of 
the mass-temperature relation can be used to constrain 
them, once the cosmological model is known from alter¬ 
native measurements (for instance, the analysis of the 
cosmic microwave background power-spectrum). To ex¬ 
plore this possibility we allow the parameters to vary 
and compare with a maximum likelihood analysis the 
predicted correlation lengths to the results obtained us¬ 
ing (a 1 ,a 2 ,a 3 ) = (6.03,2/3,0), which is our input pa¬ 
rameter set. We recall that the corresponding correla¬ 
tion lengths are ro = 15.9±0.4 and 22.8±0.7/i _1 Mpc for 
clusters with 2 < 0.3 and 2 > 0.3, respectively. Since the 
previous results showed that the dependence on the red¬ 


shift evolution parameter a 3 is negligible, we decide to 
exclude it from the present analysis. Figure [li] shows the 
68.3 and 95.4 per cent confidence levels for the samples 
of clusters with 2 < 0.3 (right panel), 2 > 0.3 (central 
panel) and for the combination of the two samples (right 
panel). Both at low and high redshift the normalisation 
and the slope of the mass-temperature relation appear 
to be strongly correlated, allowing large regions in the 
parameter space. This is expected: the same minimum 
mass can be obtained using different combinations of a 1 
and a 2 . When we combine the results of the two sam¬ 
ples, the contours shrink and the parameters are better 
constrained. However the degeneracy between the two 
parameters does not permit their determination with 
(68.3 per cent) uncertainties less than 20 per cent. 


6 CONCLUSIONS 

In this paper we presented predictions for the clus¬ 
tering properties of galaxy clusters detectable through 
the thermal Sunyaev-Zel’dovich effect. In particular we 
showed results for surveys which are expected from the 
future Planck satellite, which will cover the whole sky 
down to a sensitivity limit of U m in = 3 x 10 -4 aremin 2 . 
Although our estimate of the sensitivity limit is quite 
conservative and very likely clusters of smaller mass or 
at higher redshift will be included into the real cata¬ 
logue, the build up of a reliable cluster sample based 
on Planck data is not a straightforward process and 
requires large efforts, starting from the simultaneous 
detection at (at least) two distinct frequencies of the 
enhancement and the decrement of the SZ signal and 
ending to the source identification and redshift deter¬ 
mination. Despite these difficulties, it is worth stressing 
that Planck will give us a unique chance to construct an 
independent and unbiased galaxy cluster sample and to 
test cosmological scenarios. Here we showed, in fact, how 
it is possible to put independent constraints on some 
cosmological parameters and/or on physical properties 
of the intergalactic medium, on the basis of such a cat¬ 
alogue. 

Our theoretical predictions were obtained for dif¬ 
ferent cosmological scenarios using a model which ac¬ 
counts for the clustering of observable objects on our 
past light-cone and for the redshift evolution of both 
the underlying dark matter covariance function and the 
cluster bias factor. A linear treatment of redshift-space 
distortions was also included. Following an approach al¬ 
ready applied to X-ray selected clusters, we make use 
of a theoretical relation between mass and temperature 
to convert the limiting sensitivity of a catalogue into 
the corresponding minimum mass for the dark matter 
haloes hosting the clusters. Based on this relation, the 
estimates for the clustering properties allow tests of the 
general scheme for the biased formation of galaxy clus¬ 
ters. 

We found that the correlation length is an increas¬ 
ing function of the limits defining the surveys. A similar 
result was obtained in previous analyses of optical and 
X-ray selected clusters. When y m i n is fixed to 3 x 1CU 4 
aremin 2 as expected for Planck, the large number of de- 
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Figure 9. The variation in the predicted correlation length ro of galaxy clusters with z < 0.3 (left panel) and z > 0.3 (right 
panel) is shown in response to changes of the parameters Qj in the mass-temperature relation. The abscissa range is normalised 
using the parameter A defined as A = (a* — afy/ia!? — aj), where and are the lower and upper limits for the parameter 
a*. The parameter ol\ is the normalisation, ranging between 3 and 10 keV; a 2 is the exponent of the mass, ranging between 0.3 
and 0.9; and <23 is the exponent of the redshift dependence, ranging between -1 and 1. Results are presented for the AC DM 
model only. The dashed regions show the 1-a ranges for ro as obtained using the parameters from Tab. hi 


z<0.3 z>0.3 Combination 



a 1 a 1 a, 


Figure 10. Confidence contours (68.3 and 95.4 per cent confidence levels; dotted and solid lines, respectively) are displayed 
for the normalisation (ay) and slope (a 2 ) of the mass-temperature relation (see eq. |j); no redshift evolution is assumed here 
(a 3 = 0). Results are given for the ACDM model only. Different panels refer to cluster samples with z < 0.3 (left panel), z > 0.3 
(central panel) and to the combination of the two samples (right panel). Filled circles show the fiducial values of ai = 6.03 and 
ol 2 = 2/3 which are the input parameters. 


tectable objects will allow a very accurate determination 
of the correlation length, with l-cr error bars of few per 
cent, much better than the existing estimates. More¬ 
over the Planck cluster sample will extend to z <; 1, 
giving the possibility to measure the redshift evolution 


of the correlation function. We found that the model 
predictions are depending on the main cosmological pa¬ 
rameters, like the matter density parameter flom and 
the power-spectrum normalisation as, while the possi¬ 
ble effect of the presence of a cosmological constant is 
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almost negligible. However, the possibility of using fu¬ 
ture clustering data obtained by Planck to constrain 
the cosmological parameters is limited by the depen¬ 
dence of the results on the model adopted for the in¬ 
tracluster medium. For example, different assumptions 
on the baryon fraction produce variations of the corre¬ 
lation length similar to those found by changing flom- 
However, the dependence on as is strong enough for it to 
be constrained, in particular with high-redshift clusters 
(see Figure ^). 

When a set of cosmological parameters is fixed, the 
cluster correlation function depends on the properties of 
the intracluster medium. We found that the Planck data 
can be used to constrain the baryon fraction with an ac¬ 
curacy of few per cent once the value of the Hubble con¬ 
stant is known. On the contrary, our results showed that 
the relation between mass and temperature, which pa- 
rameterises the whole history of the physical processes 
inside the clusters, can only be poorly determined from 
clustering data alone. 
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